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ABSTRACT 

In a volume-limited sample of 63 ultracool dwarfs of spectral type M7-M9.5, we have obtained high- 
resolution spectroscopy with UVES at the Very Large Telescope and HIRES at Keck Observatory. 
In this second paper, we present projected rotation velocities, average magnetic field strengths, and 
chromospheric emission from the Ha line. We confirm earlier results that the mean level of normalized 
Ha luminosity decreases with lower temperature, and we find that the scatter among Ha luminosities 
is larger at lower temperature. We measure average magnetic fields between and 4kG with no 
indication for a dependence on temperature between M7 and M9.5. For a given temperature, Ha 
luminosity is related to magnetic field strength, consistent with results in earlier stars. A few very 
slowly rotating stars show very weak magnetic fields and Ha emission, all stars rotating faster than 
our detection limit show magnetic fields of at least a few hundred Gauss. In contrast to earlier-type 
stars, we observe magnetic fields weaker than lkG in stars rotating faster than ~ 3kms _1 , but we 
find no correlation between rotation and magnetic flux generation among them. We interpret this 
as a fundamental change in the dynamo mechanism; in ultracool dwarfs, magnetic field generation is 
predominantly achieved by a turbulent dynamo, while other mechanisms can operate more efficiently 
at earlier spectral type. 

Subject headings: stars: low-mass, brown dwarfs - stars: magnetic fields 



1. INTRODUCTION 

Ultracool dwarfs are objects of spectral types M7- 
M9.5. They can be stars in the mass range 0.08- 
0.09 M on the Main Sequence, young brown dwarfs, 
or even very young objects of planetary mass. A large 
fraction of ult racool dwarfs sho w strong chromospheric 
emission (e.g., W est et ail [2004') probably due to mag- 
netic fields that are sustained by a stellar dynamo 
Reiners fe Bas ri 2007). At masses lower than 0.3 M Q 



Mohantv & Basr| (f2003h . and iReinersI ((20071 ). While 
Reid et al.l (120021) report no relation between rotation 



..Burrows et al. 19971). ul tracool dwarfs are fully convec- 
tive ( Baraffe et al.l I1998T ) and cannot harbor a dynamo 



that requires a tachocline, as is believed to be the case 
at least for the c yclic part of the solar dynamo (e.g., 
lOssendriivenl2003D . 

The relation between rotation and activity in solar- 
type st ars and early-type M dw arfs was invest igated, 
e.g., bv lPizzolato et all (|2003f) and|ReinersJ (|2007f) . They 
find that in slowly rotating stars, i.e., stars with Ro = 
-P/tcoiw 0.1, with P the rotation period and r conv the 
convective overturn time, activity is stronger at more 
rapid rotation. Above this threshold, activity becomes 
saturated and does not grow as the stars rotate more 
rapidly. For early- and mid- type M stars (M2-M6), 
iReiners et al.l (|2009a[) showed that magnetic fields also 
saturate at this threshold. 

Among ultracool dwarfs, the exist ence of a r o tation - 
activity connection was probed by Re id et al.l (2002), 



Electronic address: Ansgar.Reiners@phys.uni-goettingen.de 
Electronic address: basri@berkeley.edu 
Emmy Noether Fellow 



and activitv. lMohantv fc Basri (2003) showed a rotation- 
activity connection down to spectral type M8 that breaks 
down at spectral type M9. They also report an in- 
crease of the saturation velocity at mid-M spectral types. 
IReiners fc Basr i (2007) report the first direct measure- 
ments of photospheric magnetic fields from Zeeman line 
broadening, but their sample is too small for a de- 
tailed investigation of the relation between magnetic field 
strength, temperature, and rotation. 

The temperature range of ultracool dwarfs is partic- 
ularly interesting for a number of reasons. Depending 
on age, stars and brown dwarfs share the same tempera- 
ture range so that young ultracool dwarfs can be brown 
dwarfs while these have evolved to later spectral type at 
higher ages. Thus, old ultracool dwarfs are always stars. 
It is also this temperature range where the atmospheric 
conditions affecting the high-energy processes of stellar 
activity are rapidly changing. At lower temperature, the 
atmosphere becomes more and more neutral, which can 
make a coupling between magnetic field lines and the 
atmosphere more difficult (Mover fc Mever-Hofmeisterl 
[ilMlMohantv et al.ll2002T) . 

In order to investigate the physics of ultracool dwarfs, 
we observed a large sample of objects with spectral types 
between M7 and M9.5. We have introd uced the sample 
in a first paper ((Reiners fc Basril 120091 Paper I in the 
following) , where we searched for Li absorption lines and 
derived space motion and the kinematic velocity of the 
sample. In this paper, we focus on rotation, activity, and 
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magnetic field measurements. 

2. SAMPLE SELECTION AND OBSERVATIONS 
2.1. The sample 

We constructed a volume-limited sample (d < 20 pc) of 
known M stars of spectral type M7-M9.5. The targets 
are taken from se veral catalogues and discoveries from 
the DENIS survey (IpeTfosse et al.ll200lHCrifo et alj |2005t 
iPhan-Bao fc Besseill2006t IPhan-Bao et Tal.1120061) Mid the 
2MASS survey (ICruz et al.ll2007D . 

From the two surveys, we constructed a joint sample 
of objects within 20 pc (for a more detailed description, 
see paper I) . Both surveys together cover more than 50% 
of the sky and can be considered as almost complete in 
the spectral range M7-M9.5. The full "20pc-2MASS- 
DENIS" survey contains 63 objects, 4 of which are found 
in both the 2MASS and the DENIS surveys. 

2.2. Observations 

Observations for our 63 sample targets were collected 
using the HIRES spectrograph at Keck observatory for 
targets in the northern hemisphere, and using UVES 
at the Very Large Telescope at Paranal observatory for 
targets in the southern hemisphere (PIDs 080.D-0140 
and 081.D-0190). HIRES data were taken with a 1.15" 
slit providing a resolving power of R « 31,000. The 
three HIRES CCDs cover the spectral range from 570 
to 1000 nm in one exposure. UVES observations were 
carried out using a setup centered at 830 nm cover- 
ing the spectral range 6400-10, 200nm with a 1.2" slit 
(R ss 32, 000). All data provide the Ha line as well as the 
absorption bands of molecular FeH around 1 /im that are 
particularly useful for the deter mination of rotation an d 
magnetic fields in M dwarfs (see lReiners fc Basril l2006'). 

Data reduction followed standard procedures including 
bias subtraction, 2D flat-fielding, and wavelength cali- 
bration using ThAr frames. HIRES data was reduced 
using the MIDAS echelle environment. For the UVES 
frames, we used the ESO CPL pipeline, version 4.3.0. 
We refer to Paper I for more details. 

3. MEASUREMENTS 

3.1. Ha Emission 

Our measure ment of Ha emission f ollows the proce- 
dure outlined in IReiners fc Basri (2008]). To measure the 
equivalent width in the Ha line against the continuum, 
we normalize the line at two footpoints blue- and redward 
of Ha. The footpoints are the median values at 6545 - 
6555 A on the left hand side, and at 6572 - 6580 A on the 
right hand side of the Ha line. None of the emission lines 
found in our targets extends into the region used for nor- 
malization. The Ha equivalent width is then measured 
by integrating the flux from 6555 to 6572 A. Equivalent 
widths of Ha emission lines are reported together with 
results on rotation and magnetic fields in Tabled The 
typical uncertainty is about a fifth of an Angstrom. We 
convert the measured Ha equivalent width into Ha flux, 
-Pkcn by multiplying it with the flux per unit equivalent 
width from the continuum flux in synthetic spectra. We 
divide the Ha flux by the bolometric flux, Fboi = erT 4 , 
and use F^ a / F\, a \ = Ln a /Lb \ to derive t he normalized 
Ha luminosity (see iReiners fc B asri 2008]). To find the 



temperature of our targets, we used the spectral types 
and calculated T e g acco rding to the conversion given in 
iGolimowski et all (|2004f ). 

3.2. Rotation and Average Magnetic Field 

The analysis of our spectra follow s the strategy laid 
out in IReiners & Basril (|2006t l2007f) and IReiners et "all 
(|2009aD . To measure the projected rotation velocity 
v sini and the magnetic field average over the visible 
surface, Bf (with / < 1), we utilize the absorption band 
of molecular FeH close to 1 ^tm. We compare our data 
to spectra of the slowly rotating M-stars GJ 1002 (M5.5) 
and G1873 (M3.5). The average magnetic flux of G1873 
was measured to be Bf = 3.9kG using an atomic Fel 
line (jJohns-Krull fc Valentil I2000D . In order to match 
the absorption strength of the target spectra, the inten- 
sity of the FeH absorption lines in the two comparison 
spectra is adjusted according to an optical-depth scal- 
ing (see IReiners fc Basril I2006D ; the FeH intensity is a 
free parameter in our fit. Note that in the magnetic 
field measurement of our template star, G1873, Bf is 
the weighted sum of several magnetic components used 
for the fit (££?/). In our data, we cannot follow such 
an approach because (1) the spectral resolution of our 
data does not allow the differentiation of individual mag- 
netic components, and (2) no information on the mag- 
netic splitting pattern of FeH lines is used. Thus, with 
our method we can only scale the product Bf from the 
spectrum of G1873. For a detail ed discussion of the sys - 
tematic uncertainties, we refer to lReiners fc Basril (|2006f) . 
The typical uncertainty in v sini is ~2kms _1 for the 
slow rotators and ~ 10 % in the case of rapid rotation. 
The uncertainty in Bf is usually several hundred Gauss. 

For the determination of the average magnetic field, 
Bf, we concentrate on relatively small wavelength re- 
gions that contain absorption lines particularly useful for 
this purpose, i.e., regions that contain some magnetically 
sensitive as well as magnetically insensitive lines. We de- 
termine Bf of our target stars by comparison of the spec- 
tral regions at 9946.0-9956.0 A and 9972.0-9981.0 A. We 
show an example of the data and the fit quality together 
with a x 2 map in v sinz and Bf in Fig.[T] In the upper 
two panels, we plot the data and two extreme cases (a) 
without any magnetic field, and (b) strong magnetic flux 
with Bf ss 4kG. Template spectra are artificially broad- 
ened and scaled to match the absorption depth of the 
FeH lines in our sample targets. The lower panel of FigfT] 
shows a % 2 -map, i.e., the goodness of fit, % 2 , as a func- 
tion of projected rotational velocity, v sini, and average 
magnetic field, Bf. We mark the formal 3tr-region for 
two free paramete rs, x 2 = x 2 lin + 11.8, with a white line 
(|Press et al.lll992t) . To determine our 3cr-uncertaintics 
in v sinz and Bf, we are marginalizing over the other 
component so that we can chose the range \ 2 — xLin + 9 
because the other parameter was varying freely. The val- 
ues of the reduced x 2 , X^i are on the order of unity except 
for the spectroscopic binary, 2MASS J0435161 - 160657, 
for which no rotation and magnetic flux analysis could 
be carried out. In one case, 2MASS J0320596+185423 
(LP412-31), Bf is probably larger than the maximum 
limit of our method (3.9 kG). Fit quality is therefore not 
optimal in this case, but the differences are likely to be 
explained by the effect of a stronger field (see also Fig. 8 
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in iReiners fe Basri|[2007f) . 

The results of our analysis for all stars are given in 
Table[3l Some of the stars were already analyzed in 
IReiners fc: Basril (j2007l ) but have been consistently rean- 
alyzed following the outlined strategy together with the 
new targets. Values of v sini and Bf are best fits from 
our x 2 analysis with formal 3tr-uncertainties. The uncer- 
tainties in Bf are usually on the order of a few hundred 
Gauss. The typical uncertainty in v sini is on the order 
of 10% with a minimum uncertainty of 2kms . Our 
detection limit due to lim ited spectral resolving power is 
v sini m i n = 3kms _1 (see IReiners fc Basri 20Q3). 

4. RESULTS 

Projected rotational velocity, v sini, normalized Ha 
luminosity, logLna/ibob an d average magnetic field, 
Bf, are plotted as a function of spectral type in Figs. [21 
[31 and and are discussed in the following. Li brown 
dwarfs as well as stars with a space velocity component 
V < — 30kms _1 are indicated (see Paper I). The lat- 
ter sample probably consists of predominantly relatively 
old objects, altho ugh this is not necessarily true for each 
individual object (|Wielenlll977l ). 

4.1. Rotation 

We find rotation velocities up to v sini » 70kms~ 1 in 
our sample. The Li brown dwarfs exhibit relatively high 
rotation velocities (with a mean of v sini = 33kms _1 ) 
while rotation in the old sample is very low (mean 
v sini = 8kms _1 ). This is a clear indication that rota- 
tional braking slows down stars of spectral type M7-M9.5 
as they age. Nevertheless, the stars of the old sample also 
exhibit significant rotation, which stands in contrast to 
old early-M type dwarfs. The longer timescales of ro- 
tational braking appear very clear in this sample (see 
IWest et al.ll200a IReiners fc Basrill2008f) . 

In the left panel of Fig. [4] we show the fraction of 
"rapid" rotators, i.e. stars with v sini > 5 km . The 
fraction of stars in our local sample populating this ve- 
locity range is about 60 % at spectral type M7. At spec- 
tral class M9, we find no object that is rotating slower 
than Skms^ 1 . Statistically, we cannot draw a robust 
conclusion from this limited sample, i.e., we cannot re- 
ject the hypothesis that the fraction of stars rotating 
more rapidly than v sini > 5km -1 is constant between 
M7 and M9.5. However, the trend that late-M dwarfs 
seem to be generally rotating more rapidly than early- 
and mid-M dwarfs is supporte d by other samples ex- 
tending into the L- dwa rf regime (jMohantv fc B asri 2003; 
IReiners fc Basril [20081 ). and seems to be rather robust. 
Such a trend is an important constraint for surveys 
searching for radial velocity variations with the goal to 
find planets around low-mass stars, because the achiev- 
able radial velocity accuracy strongly depends on rota- 
tional broadening. 

4.2. Ha-Activity 

While rotation of ultra-cool dwarfs is braked with time, 
no age effect is visible in chromospheric Ha activity 
(Fig. [3]). It is known that the mean level of normalized 
Ha activity is diminishing from mid-M to late-M and 
L spectral classes. In our sample, we observe a mean 
level of logLHa/^boi = —4.2 at spectral type M7, and of 



logiHa/^boi = —4.8 at spectral type M9. Interestingly, 
we see no clear dependence of normalized Ha on age; we 
find the few Li brown dwarfs to be on the average level 
of Ha luminosity at spectral type M8 and at the lower 
activity range at M9. Note that the high Ha emission in 
the Li brown dwarf at M7.5 is probably due to accretion. 

In the right panel of Fig.[U we plot the fraction of 
active stars among our sample (with la- uncertainties). 
Among the 63 objects, we found only three objects with- 
out any sign of Ha emission. All other 60 objects exhibit 
Ha in emission. The activity fraction in our sample is 
close to 100% at spectral types M7-M8.5 (98^ %, only 
one out of 41 shows no Ha emission). On the other hand, 
we find a fair fraction of stars with no Ha emission at 
spectral class M9-M9.5 (83^5 %). It is important to re- 
alize that the non-detection of Ha in our data does not 
exclude the presence of activity altogether; at some low 
level, all stars might show Ha emission. Again, from our 
sample it is not possible not draw robust statistical con- 
clusions about a dependence of activity on spectral type, 
but a t urnover of the ac tivity fraction at spectral type 
M9 (cp. IWest et alJl2008ft seems to be very plausible. 

4.3. Average Magnetic Field 

In the magnetically sensitive absorption lines of molec- 
ular FeH, we can measure average magnetic fields in rel- 
atively rapidly rotating M dwarfs. We report magnetic 
field measurements for all stars of our sample with a pro- 
jected rotation velocity of v sini < 20kms _1 . Above 
that rotation rate, spectral lines are too broad for a reli- 
able detection of Zeeman broadening. In all 46 stars that 
fulfill this criterion, our 3 a limits allow a determination 
of Bf within an uncertainty less or equal to 1 kG, with 
a much smaller uncertainty in the slow rotators. 

We detect non-zero magnetic fields in 41 of our 46 tar- 
gets, i.e., only in 5 (11 %) of our objects, the average mag- 
netic field is on the order of only a few hundred Gauss or 
less. Our measurements of Bf are plotted as a function of 
spectral type in Fig. [5] At all spectral types, the scatter 
among average magnetic field values is quite large. We 
find both, weak and strong fields in M7, M8, and M9 ob- 
jects. The mean average magnetic field is Bf ss 1600 G 
with an rms of 900 G. The mean of our 3<r uncertain- 
ties is ~ 430 G implying that the scatter in Bf is much 
larger than the formal statistical uncertainties. There is 
some evidence for a lack of very weak fields among the 
M8.5-M9.5 stars; all 15 stars in this spectral range show 
average fields on the order of 1 kG or stronger. Inter- 
estingly, the only object with a very weak field is the Li 
brown dwarf 2MASS J0443376+000205. On the other 
hand, six out of 30 stars with spectral types M7.0-M8.0 
have average magnetic fields weaker than 900 G. Perform- 
ing a K-S test on the two distributions among the stars 
only, we find that the hypothesis of a single underlying 
distribution for the average magnetic fields in the M7- 
M8 and M8.5-M9.5 spectral bins can only be rejected at 
the la-level (and the level is below la if we include the 
b rown dwarf). 

IReiners et al.l (|2009bf ) found that accreting young 
brown dwarfs seem to have very weak fields, much weaker 
than accreting stars at the same age (but earlier spectral 
type). There is mounting evidence that it is difficult to 
produce magnetic fields of kilo-Gauss strength at young 
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ages in brown dwarfs. Whether this applies to brown 
dwarfs in general, or is due to the young age of very-low 
mass objects, is not yet clear and can only be decided 
from magnetic field observations of older brown dwarfs, 
which are currently not available. 

Other than this, we see no evidence for a dependence 
between spectral type and average magnetic field, in par- 
ticular no rule for the production of very strong fields. 
Furthermore, there is no clear evolution of magnetic field 
strength with age. 

4.4. The Connection between Rotation, Activity, and 
Magnetic Flux 

Solar-type stars as well as early-M stars show a 
well-documented relation between rotation and mag- 



netic activity (e.g., Pizzolato et aTl [2003: Re iners! [20071 : 
iReiners et al.ll2009a ). This rot at ion- activity relation is 
generally explained by stellar activity being proportional 
to magnetic flux generation. The latter is ruled by the 
Rossby number, Ro = P/t coi1v , the ratio between the 
rotation period and the convective overturn time. For 
Ro > 0.1, magnetic flux generation and activity are pro- 
portional to the stellar rotation rate. Below that thresh- 
old, activity and magnetic field generation are saturated 
and do not grow with decreasing Ro. 

Mainly because of the difference in radius between sun- 
like and M-type stars, the threshold surface rotation ve- 
locity is a few ten kms -1 in sun-like stars, but only on 
the order of lkms -1 in mid- and late-M dwarfs. Thus, 
essentially all M dwarfs with detectable rotation (typi- 
cally v sini > 3kms _1 for m ost spectrogra phs) are in 
the saturated regime (compare iReiner si 12001 . The rela- 
tion between Rossby number and activity has been ex- 
pected to break down in fully convective stars, because 
the main motivation for the Rossby dependence comes 
from the assumption of an afi-dynamo that is proba- 
bly only operating in the prese nce of a radiative core 
and not in fully convective stars (jDurnev et~aLl[T993T ). It 
is therefore assumed that in mid- to late-M dwarfs, the 
rotation-activity relation is no longer valid or changes 
substantially. 

Three high-resolution surveys investigated the 
rotation-activity relation beyond the threshold to com- 
plete conv ection, i.e., in obj ects of spectral type mid-M 
and later. iReid et all (2002) find no connection between 
rotation and activ i ty in their sample of M7-M9.5 stars. 
iMohantv fe Basrl (|2003l ) show that a rotation-activity 
connection exists down to spectral type M8.5, and that 
the connection is weaker at M9 until it breaks down 
alt ogether among th e L dw arfs. This result is supported 
by IReiners fcF jasri (2003) who show that a T-shaped 
rotation-activity relation exists in stars earlier than 
spectr al type M9. The ir samp le com bine the objects 
from Mohantv & Basri (2003) and IReiners fc Basril 
(2008) showing that no relation between rotation 
and activity appears a t spectral types M9 and later. 
IMohantv fc Basrl (|2003l ) add that there appears to be a 
change in the magnitude of the threshold velocity above 
which Lna/iboi saturates. While this velocity is about 
3kms _1 at early-M, it appears to be roughly lOkrns -1 
at spectral types late-M according to that work. 

In our new volume- limited sample of M7-M9.5 dwarfs, 
we can now investigate the transition from the regime 



TABLE 1 

Correlations between rotation, Ha- activity, and 

THE AVERAGE MAGNETIC FIELD Bf. 



Relation 



M7 



MS 



M9 a 



Correlation coefficients 



logL Ha /L bo i vs. v sini 0.63 0.32 -0.34 (-0.13) 
Bf vs. v sini 0.32 0.69 -0.09 (-0.09) 

logi HQ /iboi vs. Bf 0.77 0.66 0.61 (0.61) 



Spearman rank test b 



logLHo/^bol vs. v sini 


0.62 


0.17 


-0.13 (0.01) 




0.006 


0.52 


0.62 (0.97) 


Bf vs. v sini 


0.32 


0.41 


0.06 (0.06) 




0.21 


0.16 


0.84 (0.84) 


logLW^bol vs. Bf 


0.58 


0.56 


0.56 (0.56) 




0.01 


0.05 


0.02 (0.02) 



a Numbers in parentheses are excluding M9 objects faster 
than v sini = lSkms -1 (see text). 

b The upper row always gives the rank-order correlation 
coefficient, r B . The lower row gives the probability of a 
non-correlation. 

of an intact rotation-activity relation towards the regime 
where this relation no longer applies. Furthermore, we 
can look for a change in the threshold velocity above 
which activity saturates. 

We divide our sample of 62 dwarfs with v sini mea- 
surements into three subsamples of consecutive spectral 
types. To obtain samples of approximately same size, 
we divide the samples into objects of spectral type M7- 
M7.5 (23 stars), M8.0 (19 stars), and M8.5-M9.5 (20 
stars). For simplicity, we refer to these subsamples as 
M7, M8, and M9, respectively. Normalized Ha activ- 
ity and average magnetic field for each subsample are 
plotted as a function of projected rotation velocity in 
Fig-El As a quantitative measure for the connection be- 
tween rotation, Ha-activity, and average magnetic field, 
we calculate the linear cor relation coefficien ts and the 
Spearman rank correlation ([Press et al.l ll992. , Chapter 
14.6) between these parameters for the three subsamples. 
The coefficients are given in Table[T] While the correla- 
tion coefficients cannot present a statistically meaningful 
measure of the significance of a correlation (particularly 
not in the case of our small sample size), the Spearman 
rank test provides a statistically more robust indicator. 
For the Spearman rank tests, we give the rank-order cor- 
relation coefficient, r s , in the first column, and the as- 
sociated probability of a non-correlation in the second 
column of Table[TJ 

4.4.1. Rotation and Ha- Activity 

Evidence exists from different earlier works that the 
rotation-activity relation undergoes a change around 
late-M spectral types. The new sample presented here 
allows for the first time to study this relation in detail in 
the narrow range M7-M9.5 and in a substantial number 
of objects that occupy a wide range in rotation velocity. 

From the upper panel of Fig. [SI the first thing we note 
is that there is evidence for an effect of supersaturation. 
In the M7 sample, objects with v sini > 30kms _1 have 
smaller values of normalized Ha luminosity than most 
of the slower rotators. At spectral type M9, the fastest 
rotators also show weaker Ha luminosity, but the satura- 
tion threshold may even be as low as v sini « 18kms _1 . 
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In the following analysis of Ha activity and rotation, 
we choose to exclude the most rapid rotators (v sin i > 
30kms _1 ). We will come back to the supersaturation 
effect in Sect. 14.51 Unfortunately this threshold velocity 
is faster than can be accommodated by our measurement 
technique for average magnetic fields. We therefore can- 
not measure Bf for the super-saturated stars. 

The linear correlation coefficients between Ha- 
luminosity and v sinz are 0.63, 0.32, and —0.34 for the 
M7, M8, and M9 subsamples, respectively (see TableQ]). 
If we exclude M9 objects faster than v sinz = 18kms _1 , 
the coefficient is —0.13. The results of the Spearman 
rank tests clearly show that clear evidence exists for a 
correlation between logLHa/^boi and v sinz in the M7 
subsamplc. On the other hand, only weak indication 
exists for such a relation at M8, but this hypothesis can- 
not be rejected either. At M9, a correlation between 
rotation and Ha luminosity is very improbable in stars 
slower than v smi — 18 kms , and if we include the 
faster rotators, rotation and activity may even be anti- 
correlated, but evidence for the latter is also weak. It is 
not surprising that including the low-activity rapid rota- 
tors at spectral type M9 leads to a weak anti-correlation. 
At this point, we cannot decide whether the drop in Ha 
luminosity around v sinz — 18kms _1 is due to a super- 
saturation occuring at this threshold velocity, or whether 
activity gradually weakens with faster rotation. 

A well-known behaviour of the normalized Ha lumi- 
nosity is that it becomes smaller with later spectral 
types. Fig. [7] shows the distribution of logLna/^boi in 
the three subsamples of objects that are rotating slower 
than v sin i = 30 km s _1 . If we assume that the difference 
between Ha luminosities within a spectral bin are not 
due to different rotation velociteis but are purely statis- 
tical, we might expect that Ln a /L^ i, or logl/Ha/^boh 
simply scatters around a mean value that depends on 
spectral type. In this case, the distribution of the nor- 
malized Ha luminosity could be described by a Gaus- 
sian fit to the data. We have tried a Gaussian fit to 
the distributions of logLua/iboi in the three subsam- 
ples; the results are overplotted in Fig. [71 Note that 
the choice of logLna/iboi (instead of Ln a /L^ i) has no 
physical motivation. In all cases, the values can consis- 
tently be described by a Gaussian distribution, which 
shows that the scatter within a spectral bin may be 
purely random. Only two outliers have very low nor- 
malized activity levels exceeding the Gaussian distribu- 
tions (one each in the M7 and M8 subsamples). Both 
are very slow rotators that probably lie on the un- 
saturated part of a rotation-activity relation similar to 
early-M dwarfs, i.e., a relation with a saturation veloc- 
ity of a few kms -1 . Two objects of the M7 subsam- 
ple show extraordinary strong Ha emission, one of them, 
2MASS J07522939+161215, was probably observed dur- 
ing a flare. The object with the strongest Ha emission 
in our sample, 2MASS 0041353-562112, is a young M8.0 
brown dwarf with detected Li. The width of the Ha line 
at 10% of the total intensity is 200 kms -1 , i.e., much of 
the observed Ha emi ssion is probably due to accretion 
(IWhite fc Basri 2003). This object is studied separately 
in lReinersI (|2009ft . 

We conclude that the scatter in Ha activity is mostly 
due to an intrinsic, purely statistical scatter that is inde- 



TABLE 2 

Mean and width of the normalized Ho 
luminosity distribution in stars with 



v smi 


< 30 KM s 


-i 


M7 


M8 


M9 a 


Mean -4.16 


-4.30 


-4.75 (-4.71) 


1 cr width 0.21 


0.36 


0.46 ( 0.42) 



a Numbers in parentheses are excluding M9 
objects faster than v sini = 18 kms -1 (sec 
text) . 



pendent of rotation and can be described by a Gaussian 
distribution. Outliers can probably be explained by very 
slow rotation, flaring, or accretion. We summarize the 
parameters of the Gaussian fits to the logLHa/iboi dis- 
tribution in the three subsamples in Tabled and we over- 
plot the 1 a range as grey shaded areas in the top panel of 
Fig.[6l While the mean level of Ha activity is decreasing 
towards later spectral type, the width of the distribution 
is growing, covering a larger area in the logLHa/^bol - 
v sinz diagram. Thus, objects like the two M8.0 stars 
at v smi rj 10 kms -1 and logLfja/iboi ~ —5 may oc- 
cupy the low end of the scatter among saturated objects 
instead of an unsaturated branch extended to higher ro- 
tation. In other words, we see no evidence for a higher 
saturation velocity at very low masses but larger scat- 
ter and lower mean values of saturated activity among 
ultracool dwarfs. 

The two main conclusions of this part are that the 
rotation-activity relation does not remain positive above 
v smi =20-30 kms -1 in any spectral type, and that 
there is not a tight correlation between velocity and ac- 
tivity in the velocity range we can resolve. 

4.4.2. Rotation and Magnetic Field 

Average magnetic fields, Bf, of stars with v smi < 
20 kms -1 are plotted as a function oft; sinz in the lower 
panel of Fig. [HI The linear correlation coefficients be- 
tween Bf and v sinz are 0.32, 0.69, and -0.09 for the 
M7, M8, and M9 subsamples, respectively. The Spear- 
man rank test shows that there is weak (la) evidence 
for a correlation between the average magnetic field and 
rotation among the M7 and M8 subsamples, and there is 
probably no correlation at M9. The strongest fields ob- 
served in our sample, average magnetic fields on the order 
of 4 kG, are observed only around v sin z of 10 km s -1 . At 
least among the M7 and M8 subsamples, very low mag- 
netic fields below our detection threshold are only seen 
in slowly rotating stars (the one M9 object with no de- 
tectable field is a young brown dwarf, see above). If we 
take away the highest and the lowest magnetic field mea- 
surement in each subsamplc, the correlation coefficients 
are 0.24, 0.12, and 0.11 (r s = 0.25, 0.12, and 0.21 with 
probabilities 0.37, 0.73, and 0.46) for the samples M7, 
M8, and M9, respectively. This means that the correla- 
tions between Bf and v sinz, that may exist in the M7 
and M8 subsamples, are mainly driven by two extreme 
values in each subsample. Other than that, the scatter in 
Bf is very large and no clear correlation is found between 
magnetic field strength and rotation. This indicates that 
at M7 and M8, very low magnetic fields still can only 
exist at very slow rotation, and that the strongest fields 
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require some rotation, but that magnetic field genera- 
tion is not as strongly tied to rotation as in earlier-type 
M dwarfs. 

iReiners et al.l (|2009a| ) showed the relation between av- 
erage magnetic field strength and rotation for a (non 
volume-limited) sample of 24 M-stars with spectral types 
between M2 and M6. In their Fig. 5, a clear relation be- 
tween v sini and Bf is visible: All stars rotating at a 
detectable rate, i.e., faster than « 3kms _1 , show av- 
erage magnetic fields stronger than k, 1500 G, and nine 
stars rotating slower than the detection threshold have 
fields below 1000 G. The conclusion for the early- to mid- 
M stars is that a rotation-magnetic field relation exists 
in the sense that all stars rotating at a detectable rate 
produce strong magnetic fields, and that no rapid rota- 
tors with weak fields are found. In the lower panel of 
Fig-El we have highlighted the region of detectable rota- 
tion and relatively weak magnetic fields (hatched area); 
this region is not populated in st ars of spectral type <M7 
(note that IReiners et alj l2009al show that this correla- 
tion probably holds for earlier stars as well) . At spectral 
type M7, there are already quite a few stars populating 
this area. At M8, field strengths in stars rotating between 
v sin? = 5 and lOkrns -1 are perhaps even weaker. Ob- 
viously, there is a marked change in the magnetic field 
strengths produced in rotating stars. This change occurs 
at spectral type M6/M7. 

In ultracool dwar fs, rotational braking is weaker than 
in earlier type stars ([Reiners fc Basril20 08). and our sam- 
ple contains only very few objects rotating at a very low 
rate. There is some evidence that very slow rotators, 
at least in the M7 and M8 subsamples, can have very 
weak fields. At detectable rotation, however, field gen- 
eration is not as efficient as in earlier-type stars. While 
there may be a weak correlation between rotation and the 
weakest magnetic fields generated at spectral types M7 
and M8, the relation probably breaks down altogether 
at spectral type M9. The main conclusion is that the 
saturation of magnetic field generation, observed in sun- 
like stars down to mid-M spectral types, breaks down 
around spectral type M7, and the efficiency of magnetic 
field generation probably depends on parameters other 
than rotation. 

4.4.3. Ha-Activity and Magnetic Field 

As we have seen in the last two Sections, the correlation 
between average magnetic flux at M7 and M8 and rota- 
tion on one side, and between Ha luminosity and rotation 
on the other side, are rather weak. We now ask the ques- 
tion whether chromospheric emission observed in Ha is 
still caused by magnetic flux, in other words, whether Ha 
luminosity still correlates with magnetic fields in ultra- 
cool dwarfs. The linear correlation coefficients between 
Ha luminosity and Bf for M7, M8, and M9 are 0.77, 
0.66, and 0.61, respectively. The Spearman rank test ex- 
cludes the null hypothesis of no correlation at the 2a level 
at all spectral types. This means that for the entire M 
dwarf range, there is evidence for a correlation between 
Ha luminosity and average magnetic field. 

We show in Fig. [8] the relation between Ha luminosity 
and average magnetic field. The three subsamples M7, 
M8, and M9 are shown with different colors and symbols, 
linear fits for the three subsamples are overplotted. In all 
three subsamples, there is a linear trend between Ha and 



magnetic fields, higher magnetic fields produce more Ha 
emission. The correlations between magnetic field and 
Ha luminosity have similar strength in all three subsam- 
ples, but they are offset with respect to each other so 
that cooler objects in general produce less Ha emission 
at the same magnetic field strength. This is probably 
explained by the growing neutrality of the atmosphere 
which makes the coupling between the magnetic fields 
and the atmosphere less efficient at lower temperatur e 
(|Mever fc Mever-Hofmeisterll 19991 : iMohantv et~ail2 002) . 

Fig.p can be compared to Fig. 8 in IReiners et alJ 
(2009a), where a similar plot is shown for spectral types 
M2-M6. At earlier spectral types, there is evidence for 
a linear relation between magnetic field strength and 
log L Ha /Lha\ f° r field strengths in the range 0-2 kG. 
Above 2 kG, Ha emission seems to be saturated and does 
not grow with higher fields. From Fig. [8) it is difficult to 
assess whether a similar saturation may occur in ultra- 
cool dwarfs because the scatter even among the individ- 
ual subsamples is substantial, and only very few objects 
occupy the region of fields higher than 3kG. We con- 
clude that Ha emission still depends on magnetic field 
strength in a way that is similar to earlier stars. The 
scatter among Ha emission at constant field strength is 
substantial, it is larger than at earlier spectral types, 
which may be due to differences in atmospheric temper- 
ature and ionization. It is currently unclear whether a 
saturation of Ha emission occurs around field strengths 
of 2 kG like in earlier stars; more observations of stars 
with very strong fields (> 3 kG) are required to answer 
that question. 

4.5. Evidence for super saturation? 

Late-M dwarfs have small radii, relatively short ro- 
tation periods, and probably large convective overturn 
times. While this renders a study of activity difficult at 
Rossby numb ers on the order of 1 (w eakly active or inac- 
tive stars, see iPizzolato et aLll2003l) . it opens the oppor- 
tunity to investigate the other extreme of the rotation- 
activity connection, i.e., to probe activity at very low 
Rossby numbers where so-called "supersaturation" may 
exist. This effect ma y be a result of coronal stripping due 
to centrifugal forces (Jjardinc 2004) or it could be due to 
some effect within the dynamo itself. 

To calculate the Rossby number, Ro = P/r conv , we 
require knowledge of the rotation period P and the con- 
vective overturn time t c . Unfortunately, no detailed cal- 
culations of the convective overturn time for all M spec- 
tral types are available, and it is even unclear where in 
the convection zone an effective overturn time should be 
defined. Furthermore, M dwarf rotation periods are diffi- 
c ult to measure. A sizeable number of periods is reported 
in lKiraga fc Stepier] (|2007h . but still, for most M dwarfs, 
no rotation periods are available. On the other hand, we 
have a large sample of v sini measurements, from which 
we now estimate the rotation period to investigate chro- 
mospheric activity at very low Rossby numbers. 

In order to convert the projected rotational velocity 
v sin i to rotation period for our latc-M sample, we first 
estimate the radii of the stars. Note that this only pro- 
vides Ro/ sini, which is an upper limit of Ro. To deter- 
mine th^adhis ! _we emplo y the mass-luminosity relation 
from iDelfosse et ail (I2000J) and the ma s s-radi us relation 
at an age of 5 Gyrs from iBaraffe et all ([1998) . To com- 
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pute t he masses, we use J-magnitudes from iCutri et al.l 
(|2003f ) and distances given in Paper I. We augment our 
sample by inclu ding the M0-M9 sta rs with v sin i mea- 
surements from IDelfosse et"afl (11998D . iMohantv fc Basril 
(2003), and iReiners fe Basril (|2007h . For them, we es- 
timate the radius from a spec tral type-ra dius relation 
based on the values reported in lCoxl (|2000l ). and we use 
R = 0.1 R.0 for spectral type M8 and later. 

We calculate the convective overturn time as a func- 
tion of mass. We estimate this function from a fi t to th e 
empirical results presented in iKiraga fc Stepienl (|2007f ). 
but wc limit the convective overturn time to a maximum 
of T conv = 70 d consis tent with the values giv en for M 
dwarf s in |Saar (200j]), which are taken from iGillilandl 
(fj~986T) . The relation is r conv [d] = 86.9 - 94.3M/M Q for 
spectral type <M7 and T conv = 70 d for M7-M9. For the 
stars not contained in our new M7-M9 sample, the values 
for the mass are simply estimated to be the same as the 
rad ii in solar units (e.g ., 0.2 R Q corresponds to 0.2 M Q ; 
see iDemorv et al.l [2009) . Although there are certainly 
more accurate ways to estimate radius and mass, this 
rather simple approach is sufficient for our more quali- 
tative investigation of this sample, and the main uncer- 
tainty still resides in the value of the convective overturn 
time. For exam ple, if we were using the c onvective over- 
turn times from IKiraga fc Stepienl (f2007f ) for our late-M 
sample (i.e., without limiting r c to a maximum of 70 d), 
the Rossby numbers would be about 0.15 dex smaller for 
the coolest stars, which is equivalent to an uncertainty 
in the Rossby number of 30%. This difference does not 
affect our conclusions, and errors in the radius estimate 
are likely smaller than this. 

We plot the normalized Ha luminosity as a function 
of Rossby number in Fig. [9) only M dwarfs with de- 
tected Ha emission are plotted. The M dwarfs occupy 
the range between —3.0 and —1.4 in Rossby number, 
which is well within the saturation regime (see, e.g., 
iPizzolato et "all [20031. There is no obvious difference 
between our subsamples of spectral types M7-M9 and 
the other objects. However, there is a hint of a lack 
of very active objects (logiHa/^boi > —4.3) with very 
low Rossby numbers (Ro < —2.6). None of the ob- 
jects in our late-M sample occupies this region while 
many objects with higher Rossby numbers are more ac- 
tive. Nine objects with Ro < — 2.6 show relatively little 
activity (log Lhq / -^boi < —4.3). This reflects the evi- 
dence for supersaturation that we saw in Fig. [6] From 
the earlier M dwarfs, however, one star, G 165-08 (M4.5, 
v sinz = 55.5 km s -1 : IDelfosse et al~lll998l ). has a Rossby 
number smaller than —2.6 and does show stronger activ- 
ity. We divided the combined sample of all M dwarfs into 
four bins of equal size in Ro (0.4 dex) and calculated the 
median and its standard error for each bin. These values 
are overplotted in Fig.[9l they show a slight trend of lower 
Ha luminosity with lower Rossby number, but the trend 
remains inconclusive because of the large scatter and the 
low number of objects particularly in the bin with the 
lowest Ro. The Spearman rank correlation yields a coef- 
ficient of 0.20, and the null hypothesis of no correlation 
can be rejected only at the lcr-level. 

M dwarfs are very variable. The highest activity levels 
are probably connected with very strong flaring, i.e., the 
most active objects may be observed during a temporal 



maximum of their activity. Although we see no clear evi- 
dence for a supersaturation effect in our sample, the lack 
of very active objects with very small Ro - at least among 
the late-M dwarfs - gives a hint towards the effect of su- 
persaturation. A plausible explanation for such an effect 
may be that coronal stripping inhibits the occurence of 
strong flares. More observations of very rapidly rotat- 
ing M dwarfs (v s'mi > 50kms _1 ) would be helpful to 
settle this issue. In order to decide whether it is the mag- 
netic field that is weaker at very low Rossby number, or 
whether it is the reduced radius of the corona, the best 
next step would be to find a way to actually measure 
the magnetic fields on the stars with very low Rossby 
numbers. 

5. SUMMARY AND CONCLUSIONS 

We investigated rotation, chromospheric emission, and 
average magnetic fields in a volume-limited sample of 63 
ultracool M dwarfs. The main results from the analysis 
are the following: (1) The mean level of Ha luminos- 
ity is decreasing with lower temperature, and the scatter 
among Ha luminosity is larger at lower temperature, but 
this scatter is not related to rotation; (2) we see a hint 
for supersaturation occuring somewhere around a Rossby 
number of —2.5, this may be caused by inhibited flaring 
due to centrifugal forces but it could also be due to dy- 
namo changes that we can't access until a method for 
measuring fields in very rapid rotators is found, and in 
any case requires confirmation from more observations; 
(3) in ultracool dwarfs of a given temperature, Ha lu- 
minosity still is a function of magnetic field strength, it 
is unclear whether the depedence shows saturation at 
high magnetic fields (as in earlier stars); (4) a few very 
slowly rotating ultracool dwarfs have very weak magnetic 
fields, but ultracool dwarfs that rotate at a detectable 
rate generate substantial fields (with one exception that 
is a young brown dwarf); (5) the saturation of magnetic 
field generation breaks down at spectral type M7, before 
M7, rapid rotation always implies a field strength on the 
order of 2 kG or higher, while at M7 and later, stars ro- 
tating as rapidly as lOkms" 1 can have magnetic fields 
on the order of 1 kG. 

The lack of saturation of magnetic fields in ultracool 
dwarfs at rotation velocities above fh 3kms _1 probably 
means that the dynamo efficiency suffers a change around 
spectral type M7. Before M7, rotation rates faster than 
P fx 2-3 d causes magnetic fields of 2kG strength or 
stronger. At M7 and later, the fields can be much weaker. 
A possible explanation for this behavior is that there oc- 
curs a change in the dynamo mechanism that predom- 
inantly generates the magnetic fields we observe. The 
reason for such a change is not at all clear. These stars 
do not harbor a tachocline, the boundary layer between 
the radiative core and the convective shell, which is be- 
lieved to play an im portant role in the solar dynamo (e.g., 
lOssendriiveri [2003V However, the tachocline is believed 
to di sappear around spectral type M3 (e.g.. iSiess et al.l 
2000), and M4-M6 stars should already be completely 
convective. So far, we see no obvious change in structure 
that could occur around spectral type M6/M7. 

Our understanding of purely turbulent dynamos is only 
at its beginning, but large progress was made during the 
last years (e.g., iDurnev et all Il993t iDobler et all [2006; 
[Browning 2008:). The observations of ultracool dwarfs 
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shown in this paper present empirical information on 
fully convective dynamos, but future investigations of 
dynamo-relevant parameters are required to put more 
constraints on dynamo modes that are realized in stars. 
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Fig. 1. — Data and fit (top panel) and \ -landscapes (bottom panel) for 2MASS 0818580+233352. In the upper panel, data are shown 
as black histograms. The three coloured lines show our fit for no magnetic field (blue line), strong magnetic flux (Bf ~ 4kG, red line), 
and the best fit, which is an interpolation of the two (green line). In the bottom panel, dark and blue colour indicates low \ 2 values, red 
and yellow colours show bad fit quality (high values of \ 2 ■ The white line shows the 3(j-level. 
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Fig. 2. — Projected rotation velocity v sini as a function of spectral type. Blue filled circles are members of the old population, red stars 
are young brown dwarfs with Li detection. Note that v sini is plotted on a log-scale, and that the old stars are predominatly occupying 
the region v sini < lOkms - ■ 
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Fig. 3. — Normalized Ha luminosity as a function of spectral type. Symbols as in Fig. [2] 
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Fig. 4. — Fraction of stars rotating at 5km/s or more (left) and fraction of active stars (right). Error bars show ler-uncertainties. 
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Fig. 5. — Average magnetic field, Bf as a function of spectral type. Symbols as in Fig. [2] Positions in spectral type are plotted with a 
small offset for clarity. 
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Fig. 6. — Normalized Ha luminosity and average magnetic field as a function of rotation velocity for different spectral types. Left, center, 
and right panels show spectral types M7-M7.5, M8.0, and M8.5-M9.5, respectively. Top panel: Normalized Ha luminosity as a function 
of spectral type. The grey shaded areas visualize the la-region around the mean for objects with v sini < 30kms — 1 (see Fig.[7jl. Bottom 
panel: Average magnetic field, Bf, as a function of rotation velocity. No measurement of Bf is possible at v s'mi > 20kms~ 1 . The 
hatched region shows the area that is not occupied by stars earlier than M7 (see Fig. 5 in Reiners ct al. 2009a). In all panels, red stars 
indicate young brown dwarfs with Li, blue filled circles show members of the old population as in Figs.|2]j5] The strong Ho; emission of the 
very active, young M7 object is probably due to accretion. 
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Fig. 7. — Histogram of normalized Ha luminosities for the three subsamples M7, M8, and M9 (top, center, and bottom row, respectively). 
Only objects with v sini < 30kms — 1 are included. Gaussian fits to the distributions are overplotted as dashed lines, the mean values and 
standard deviations of the distributions are given in Tablc[2] 




Fig. 8. — Normalized Ha luminosity as a function of average magnetic field, Bf. The three subsamples M7, M8, and M9 are shown with 
different symbols and colors (blue squares — M7; black circles — M8; red triangles — M9). A linear fit was applied to each subsample, they 
are overplotted as dashed lines. 
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Fig. 9. — Normalized H« luminosity plotted as a function of Rossb y number for our samp le (blue squares - M7; black circles — M8; red 
triangles — M9) and results taken from Dclfossc et al. (1998), Mohanty & Basri (2003), and Reiners & Basri (2007) (open circles). Green 
stars mark the median values when the sample is binned into four chunks, the error bars show the standard error of the median within 
each bin. 



